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ABSTRACT 

We perform detailed calculations of cosmic microwave background (CMB) anisotropies in a CDM-dominated 
open universe with primordial adiabatic density perturbations for a variety of reionization histories. We show 
that to a great extent, the CMB anisotropies depend only on the geometry of the Universe, which in a matter 
dominated universe is determined by O, and the optical depth to the surface of last scattering. In particular, 
the location of the primary Doppler peak depends primarily on O and is fairly insensitive to the other un- 
known parameters, such as O;,, /i. A, and the shape of the power spectrum. Therefore, measurements of CMB 
anisotropies on small scales may be used to determine O. 



Submitted to The Astrophysical Journal Letters 
kainion@guinness.ias.edu 
dns@astro.princeton.edu 
sugiyama@bkyast.berkeley.edu 



t76u^re eooioi76/Md-ojis^ 



1. INTRODUCTION 



2. REIONIZATION 



The recent discovery of anisotropies in the cosmic microwave background 
(CMB) (Smoot et al. 1992) has ushered in a new era in cosmology (for a recent 
review, see White, Scott, & Silk 1993). Although the current detections are still 
accompanied by significant uncertainties, the CMB anisotropies on all angular 
scales will soon be mapped out with great precision. 

The detailed shape of the CMB spectrum can depend quite sensitively on 
the uncertain parameters involved, and it is possible that detailed measurements 
of CMB anisotropies can be used to determine some of the uncertain parameters 
such as the Hubble constant, the baryon density, and spectral index (Bond et 
al. 1993). In this Letter, we point out that CMB anisotropies may also provide 
information on the value of Q. 

In a previous paper (Kamionkowski & Spergel 1993; hereafter referred to 
as paper I), we computed the large-angle anisotropies in a CDM-dominated 
open universe with primordial adiabatic density perturbations. We found that 
models with 0.35 < fi < 0.8 and 0.4 < /i < 0.6 (where h is the Hubble constant 
in units of 100 km sec""*^ Mpc""*^) are compatible with the observed amplitude 
of large-angle CMB fluctuations, the amplitude and spectrum of fluctuations of 
galaxy counts in the APM, CfA, and QDOT surveys, and an age of the Universe 
greater than 13 Gyr. 

The analytic treatment in Paper I is valid only for angular scales which 
probe comoving distance scales well outside the horizon at the surface of last 
scattering. On angular scales smaller than those probed by COBE, peculiar 
motions of the plasma at the surface of last scattering must be taken into 
account. An accurate treatment requires numerical solution of the Boltzmann 
equations. 

In the next Section, we discuss the reionization of the Universe by an early 
generation of star formation and argue that reasonable reionization histories 
imply an optical depth to Thompson scattering of ~ 0.5 — 1 between us and 
the redshift of recombination. In Section 3, we briefly describe our numerical 
results, and compare with our previous analytic calculations. We also show 
that the CMB spectrum depends primarily on Q and the optical depth to the 
surface of last scattering. In the final Section, we make concluding remarks. 



In the standard big-bang scenario, most of the plasma in the Universe 
recombined to form neutral hydrogen and helium at 2 ~ 1300. Yet observations 
of high redshift quasars at 2 ~ 4 imply that the intergalactic medium has been 
mostly ionized since at least that redshift (Jenkins & Ostriker 1991). This 
implies that sometime between 2 ~ 4 and z ~ 1300 the Universe was reionized 
perhaps by a burst of star formation or by early active galactic nuclei. This 
reionization at redshift zn has a significant effect on the CMB spectrum (see, e.g. 
Sugiyama, Silk, & Vittorio 1993) as it suppresses fluctuations at angular scales 
much smaller than the horizon size at that redshift, 6h ~ 2[Ci/zRY'^, and 
produces additional fluctuations both through Doppler scattering at angular 
scales comparable to 6h and at arcsecond to arcminute scales through second 
order effects (Ostriker & Vishniac 1986, Vishniac 1987). 

If we assume that the first stars are similar to those observed today, then 
we can estimate the number of ionizing photons produced per hydrogen atom 
turned into stars. Most of the ionizing radiation is produced by massive stars on 
the main sequence. A 30 Mq star burns 7 Mq of hydrogen releasing O.O5M0C'^ 
of energy. Roughly 1/4 of this energy is emitted as ionizing radiation (Kurucz 
1979). Since these massive stars account for approximately 1/4 of the initial 
mass function, star formation releases ~ 100 keV in ionizing radiation per 
processed baryon. For an O star spectrum, for a typical ionization, roughly 
13.6 eV goes to ionize hydrogen and an additional 4 eV goes into photoelectric 
heating. Thus, there are ~ 5000 ionizing photons emitted per processed baryon. 
This suggests that a minimum of 2 x 10"'* of the Universe must be in bound 
objects for the Universe to be mostly reionized. 

These ionizing photons will create a Stromgren sphere of ionized material 
around each collapsed region. The mass of this Stromgren sphere will be de- 
termined by the balance between ionization and a('^'(T), the temperature (T) 
dependent recombination rate to energy levels above the ground state (Spitzer 
1978). The number of ionizing photons needed per hydrogen atom to keep the 
gas mostly ionized at redshift z is therefore, ~ t[z)/ns{z)a''^>(T), where t[z) is 
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FIG. 1. This figure shows the fraction of the mass of the Universe in collapsed 
objects as a function of the optical depth between the present and a given red- 
shift. These calculations assumed h = 0.6, fij/i'^ = 0.015 and a scale-invariant 
spectrum of density fluctuations. The shaded region corresponds to the fraction 
of baryons in stars required to ionize the Universe. 

the Hubble time, and ne[z) is the electron density. For reasonable assumptions 
about the electron temperature, this implies that ~ 5 — 20 photons must be 
produced per hydrogen atom to keep the Universe reionized back to r ~ 1. This 
suggests that the Universe should become reionized when the fraction of mass 
in bound objects reaches a level of ~ 2 x 10"'* — 4 x 10~^. 

We can estimate the fractional mass of the Universe that has collapsed 
into bound objects of mass greater than the Jeans mass, fcoii{Mj), by using 
the Press-Schechter formalism (Press & Schechter 1974): 
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fcoiiiMj, Q, z) = -Erfc I —= '■ 
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Here, a[Mj,0) is the amplitude of mass fluctuations today at the Jeans mass 
scale and D is the linear-theory growth factor. The Jeans mass after recombi- 
nation is ~ 10^ Mq. Fig. 1 shows the fraction of mass in collapsed objects as 
a function of t{z) ~ 0.04:QihQ~^/^XeZ^/^ for different COBE-normalized cos- 
mologies. Using our estimate of required mass fraction for ionization implies 
that the Universe is reionized at a redshift of 40 - 50 creating an optical depth 
to Thompson scattering between us and large redshifts of ~ 0.5 — 1. Bardeen et 
al. (1987) argue that unreionized open models violate limits from the OVRO 
experiment on the amplitude of fluctuations on small angular size. However, 
the optical depth estimated here is sufficient to smooth out the fluctuations on 
these scales and evade this limit. As we will see in the next Section, ongoing 
experiments will enable us to constrain r for this family of models. 

3. RESULTS 

The CMB spectra are computed using the numerical techniques described 
by Sugiyama & Gouda (1992). We assume a scale-invariant Harrison-Peebles- 
Zel'dovich primordial spectrum of density perturbations. On scales larger than 
the curvature scale, the deflnition of scale invariance is somewhat ambiguous, 
but in Paper I it was shown that this ambiguity affects only the lowest multipole 
moments (Z < 9 for Q ~ 0.1 and I < 5 for Q ~ 0.3), and that the CMB spectrum 
is relatively insensitive to the exact shape of the power spectrum on scales larger 
than the curvature scale. 

In Fig. 2 we plot the COBE-normalized CMB spectrum as a function of 
multipole moment I for several values of Q and for optical depths r = (no 
reionization) and r = 1, where r ~ 0.04:CiihCi~^'^Xs{zis)^'^. Here, fij is the 
mass density of baryons, zis is the redshift of the surface of last scattering, 
and Xe is the ionization fraction. Note that the location of the flrst Doppler 
peak increases as Q is decreased, and is insensitive to the ionization history. 
Reionization decreases the amplitude of the Doppler peak by roughly a factor 
e~^'^ , and leads to another smaller peak at lower I which arises from peculiar 
motions of scatterers at low redshift in a reionized model. At large I, the 



results are consistent with Bond et al. (1993). 
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FIG. 2. The COBE-normalized CMB spectrum as a function of multipole 
moment I for several values of Q and for optical depths r = (no reionization) 
and T = 1. Here we have taken fij = 0.06 and h = 0.5. Throughout the paper, 
A = 0. 

multipole spectrum drops sharply at the angular scale subtended by the Silk 
damping scale at recombination. The location of this fall-off is sensitive to Q. 
Interestingly enough, in all these models the spectra on COBE scales, 
Z < 13, are very close to flat, especially when we consider the effect of cosmic 
variance; at Z = 10 this introduces a la error of about 30% in the prediction. 
The numerical results for the shape and amplitude of the spectrum on these 
scales agrees well with the analytical results in Paper I, especially for the case 
of no reionization. On scales I > 15, the tail end of the Doppler peak becomes 
important and a numerical calculation is needed. For fi = 1, the numerical 
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Table 1. Predictions for (AT/T)rms x 10^ for several experiments for various 
values of Q and r and the reported level of observed CMB fluctuations. The 
two values quoted for the MSAM experiment represent the r.m.s. tempera- 
ture fluctuations with and without "sources" (Cheng et al. 1993, Page 1993). 
The two values for the MAX experiment are from two different regions of the 
sky. The higher value (Gundersen et al. 1993) is reported as a detection. The 
lower value (Meinhold et al. 1993) may be fluctuation due to free-free emis- 
sion, synchrotron emission, or may represent CMB fluctuations. For the White 
Dish experiment (Tucker et al. 1993), we use the "Method I" 2-beam differ- 
encing scheme rather than the "Method II" quadrupole differencing scheme. 
The ARGO (de Bernardis et al. 1993) and SP (Schuster et al. 1993) results 
are the amplitudes of Gaussian correlation function; the r.m.s. temperature 
fluctuations will differ slightly from these numbers. 
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FIG. 3. The COBE-normalized CMB spectrum as a function of multipole 
moment I for several values of fij and h. In all models, Q = 0.3 and r = 1. 

In Table 1, we list the predictions for (AT/T)rms for the various models 
for several different CMB experiments. All of the models are computed for 
fij = 0.06 and h = 0.5. Note that the detection of CMB fluctuations at the 
level of greater than 10"^ by the MAX (Meinhold et al. 1993; Gundersen 
et al. 1993), MSAM (Cheng et al. 1993; Page 1993), SP (Schuster et al. 
1993), Saskatoon (Wollack et al. 1993), and ARGO (de Bernardis et al. 1993) 
experiments suggest that the optical depth of the Universe due to reionization 
is less than 2. This is consistent with our theoretical expectations for Gaussian 
theories. 

In Fig. 3 we plot the CMB spectrum for several values of fij and h for 
Q = 0.3, and for fixed optical depth. The similarity of the curves demonstrates 
that the primary effect of varying fij, h, zis, and x^ can be determined by their 



scaling with r. Variation of any of the parameters with fixed r has relatively 
little effect on the CMB spectrum. 

Fig. 2 and Fig. 3 show that variations in the location of the Doppler peak 
are highly insensitive to all the uncertain parameters except for the geometry 
of the Universe. Although we have used only scale-invariant spectra, the re- 
sults of Bond et al. (1993) show that the location of the Doppler peak is also 
relatively insensitive to changes in the primordial spectral index, the addition 
of a cosmological constant, or the effects of gravitational waves. If the Universe 
is flat but dominated by a cosmological constant, then the total mass density 
(non-relativistic matter plus vacuum-energy density) of the Universe is fi = 1, 
and the Doppler peak is still located at I ~ 200 (Bond et al. 1993). The ro- 
bustness of the prediction for the location of the Doppler peak as a function of 
Q should come as no surprise: This simply reflects the angular scale subtended 
by the horizon at the surface of last scattering. 

Information on the value of Q should soon be available even without fully 
mapping out the angular CMB spectrum around the Doppler peak. If fi = 1, 
then the amplitude at I ~ 200 is greater than that at I ~ 800, but if fi < 1, the 
ratio is reversed. In Figure 4, we plot the predictions of scale-invariant models 
with different ionization histories and values of Q. Figure 4 shows that all of the 
flat models tend to along a curve in this temperature ratio plot. The various 
vacuum-dominated and tilted models studied by Bond et al. (1993) also lie 
along this curve. The position along the curve depends mostly upon the optical 
depth of the Universe. Open models with Q < 0.5 lie in upper portion of the 
plot. This suggests a combination of experiments probing I ~ 200 and I ~ 700 
could potentially determine the geometry of the Universe. 

4. CONCLUSIONS 

In paper I, we showed that COBE-normalized low-fi models are consistent 
with galaxy counts in large-scale surveys and with the age of the Universe. 
Due to cosmic variance, as well as the uncertain contribution of tensor modes. 
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FIG. 4. The points in this plot represent the predicted ratios of temperature 
fluctuations for different cosmological parameters. The open symbols represent 
T = models, the star-shaped symbols represent r = 0.5 and r = 1 models, the 
filled symbols represent the r = 2 models. The triangular and three-pointed 
symbols represent fi = 0.2 models. The square and four-pointed symbols repre- 
sent fi = 0.5 models. The pentagons and five-pointed symbols represent fi = 1 
models. The solid, dashed and dotted lines connect the fi = 0.2, fi = 0.5 and 
n = 1 models. In the source-free MSAM data set, MSAM2/COBE = 1.6 ±0.5. 
However, if the sources are included MSAM2/COBE = 3.5 ± 1.0. The current 
White Dish upper limit is not yet strong enough to constrain models. 



We have also estimated when the first generation of stars may have reion- 
ized the intergalactic medium in these models. Because of the strong depen- 
dence of the number of collapse objects on redshift, the predicted epoch of 
reionization is relatively insensitive to assumption about the details of primor- 
dial stars (as long as some massive stars are part of this first stellar generation). 
In the low fi models considered here, we expect that the Universe should have 
an optical depth of ~ 0.5 — 1 between us and the epoch of recombination. This 
additional opacity to photon scattering reduces the predicted amplitude of fluc- 
tuations on small angular scales but does not shift the location of the Doppler 
peak. 

The location of the Doppler peak, I ~ 200/0-*^" reflects the angular scale 
subtended by the horizon at the surface of last scattering, so it should be 
no surprise that it is relatively insensitive to other cosmological parameters. 
Similarly, the location of the sharp drop-off in the multipole spectrum at I ~ 
9 X 10^fi~^''*(fii/i)-'^", the angle subtended by the Silk damping scale at the 
surface of last scatter, also depends on fi. Although the amplitude of the 
Doppler peak may be suppressed by reionization, it is quite likely that it will 
still be distinguishable for reasonable ionization histories. 

Current measurements are still too uncertain to discriminate between the 
various values of fi. Accurate results from the White Dish and MSAM experi- 
ments should be able to distinguish between a flat universe and a universe with 
fi ~ 0.3. Experiments with higher resolution which probe multipole moments 
100 < Z < 1000 (angular scales 0.1° < 61 < 1°) will be able to make the deter- 
mination more precise. We hope that the ideas presented here will give added 
impetus to the development of such detectors. 



discrimination between the various values of fi based on measurements of large- 
angle CMB anisotropies alone was currently impossible and likely to remain so. 
In this Letter, we have performed detailed numerical computations of the CMB 
spectrum on smaller angular scales and found that information on the value of 
Q is most likely encoded in the angular spectrum of CMB anisotropies. 
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